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Abstract  The  terminal  guidance  problem  for  an  unpowered  lifting  reentry  vehicle  against  a  sta¬ 
tionary  target  is  considered.  In  addition  to  attacking  the  target  with  high  accuracy,  the  vehicle  is 
also  expected  to  achieve  a  desired  impact  angle.  In  this  paper,  a  sliding  mode  control  (SMC)-based 
guidance  law  is  developed  to  satisfy  the  terminal  angle  constraint.  Firstly,  a  specific  sliding  mode 
function  is  designed,  and  the  terminal  requirements  can  be  achieved  by  enforcing  both  the  sliding 
mode  function  and  its  derivative  to  zero  at  the  end  of  the  flight.  Then,  a  backstepping  approach 
is  used  to  ensure  the  finite-time  reaching  phase  of  the  sliding  mode  and  the  analytic  expression 
of  the  control  effort  can  be  obtained.  The  trajectories  generated  by  this  method  only  depend  on 
the  initial  and  terminal  conditions  of  the  terminal  phase  and  the  instantaneous  states  of  the  vehicle. 
In  order  to  test  the  performance  of  the  proposed  guidance  law  in  practical  application,  numerical 
simulations  are  carried  out  by  taking  all  the  aerodynamic  parameters  into  consideration.  The  effec¬ 
tiveness  of  the  proposed  guidance  law  is  verified  by  the  simulation  results  in  various  scenarios. 

©  2014  Production  and  hosting  by  Elsevier  Ltd.  on  behalf  of  CSAA  &  BUAA. 


1.  Introduction 


For  the  case  of  the  terminal  guidance  phase,  the  guidance  laws 
are  required  to  hit  the  target  from  a  specific  direction  with  high 
precision,  as  well  as  to  achieve  a  minimal  miss-distance.  Tradi¬ 
tional  guidance  laws,  such  as  classical  proportional  navigation 
guidance  (PNG),  are  very  popular  in  practical  applications 
because  of  the  ease  of  mechanization  due  to  less  information 
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demand.  Although  such  guidance  laws  can  be  used  to  attack 
the  target  accurately,  they  are  usually  silent  on  terminal  angle 
constraint. 

Since  the  concept  of  impact  angle  guidance  is  advanced  in 
Ref.2,  great  development  has  been  achieved  on  this  issue.  A 
technique  that  has  been  commonly  employed  to  derive  such 
laws  is  using  the  PNG-based  method.  A  biased  PNG  law  is 
presented  in  Ref.3,  and  the  impact  angle  is  guaranteed  by  add¬ 
ing  a  biased  term  to  PNG.  Similarly,  in  Ref.4,  an  improved 
biased  PNG  law  which  is  more  simple  and  practical  than 
Ref.  is  proposed.  An  adaptive  guidance  is  presented  in 
Ref.5,  where  nonlinear  parameter  adaptation  laws  are  adopted 
in  PNG  to  intercept  a  ground  target.  This  guidance  approach 
is  effective  and  easy  to  implement.  Moreover,  the  trajectories 
generated  by  this  method  are  similar  to  the  optimal  trajecto¬ 
ries.  Another  PNG-based  guidance  law,  possessing  the  ability 
of  intercepting  a  stationary  target  with  various  angles  by 
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varying  the  PN  gains,  is  proposed  in  Ref.6.  Afterwards,  the 
method  in  Ref.6  is  further  extended  to  deal  with  nonstationary 
nonmaneuvering  targets  in  Ref.  .  However,  the  control  effort 
is  not  continuous  at  the  instant  of  guidance  law  switching, 
which  may  increase  the  tracking  difficulty  for  the  flight  control 
system. 

In  some  recent  work,  the  time-to-go  (tgo)  information  is 
used  to  design  the  impact  angle  guidance  laws.  A  guidance 
law  on  the  basis  of  state  dependent  Riccati  equation  (SDRE) 
is  proposed  in  Ref.  .  The  guidance  problem  turns  to  be  a  non¬ 
linear  regulator  problem  and  a  tgo  based  state-weighting 
matrix  is  used.  A  closed-form  solution  of  a  linear  quadratic 
optimal  control  problem  to  cater  for  terminal  angle  con¬ 
straint  is  developed  in  Ref.9.  In  addition,  three  kinds  of  meth¬ 
od  are  introduced  to  calculate  tgo.  However,  in  practical 
applications,  the  precise  estimation  of  tgo  is  usually  a 
challenging  work. 

On  the  other  hand,  a  computationally  efficient  technique 
which  is  called  model  predictive  static  programming  (MPSP) 
is  used  to  obtain  impact  angle  constrained  guidance  laws. 
MPSP  is  based  on  nonlinear  optimal  control  theory,  and  an 
effective  trajectory  optimization  concept  is  integrated  into 
the  guidance  laws.  For  more  details,  one  can  refer  to  Refs.10,11. 

Sliding  mode  control  (SMC),  as  an  effective  method  to 
tackle  the  uncertainties  and  disturbances,  has  drawn  a  good 
amount  of  attention  from  researchers  around  the  world.  In 
Ref.  ',  a  new  guidance  concept  named  intercept  angle  guid¬ 
ance  is  proposed.  This  guidance  method  which  is  implemented 
in  the  usage  of  sliding  mode  approach  can  be  applied  in  three 
kinds  of  engagement  scenarios,  i.e.,  head-on,  tail-chase  and 
head-pursuit.  Then,  this  guidance  concept  is  further  investi¬ 
gated  in  Ref.  3  with  taking  time-varying  acceleration  bounds 
into  consideration.  A  novel  sliding  mode-based  impact  time 
and  angle  guidance  law  is  presented  in  Ref. !4.  The  line-of-sight 
(LOS)  angle  history  is  shaped  to  satisfy  the  needs  of  impact  an¬ 
gle  as  well  as  impact  time,  then  a  backstepping  based  second- 
order  sliding  mode  control  is  derived  to  track  the  pre-designed 
LOS  angle.  The  method  presented  in  Ref.15  is  a  terminal  SMC 
(TSMC)-based  guidance  and  control  law  which  can  be  used  to 
intercept  stationary,  constant-velocity,  and  maneuvering  tar¬ 
gets  with  any  impact  angle,  even  under  the  circumstance  of 
large  initial  heading  errors.  In  spite  of  the  effectiveness  of 
the  proposed  method  from  the  constant-velocity  missile  point 
of  view,  it  is  difficult  to  estimate  its  performance  by  taking  all 
the  aerodynamic  parameters  into  consideration.  In  Ref.16,  a 
scheme  of  integrated  guidance  and  control  on  the  basis  of  an 
adaptive  SMC  algorithm  is  developed.  In  this  work,  the  non¬ 
linear  dynamics  is  taken  into  account,  and  the  robust  charac¬ 
teristics  of  the  algorithm  are  validated  in  the  case  of 
aerodynamic  parameter  perturbation;  however,  the  range  of 
available  terminal  angle  is  limited. 

In  this  paper,  an  SMC-based  guidance  law  is  developed  to 
hit  a  stationary  target  from  a  specified  direction  with  high  pre¬ 
cision.  The  nonlinear  point-mass  dynamics  in  the  pitch  plane  is 
used  to  evaluate  the  guidance  law.  As  the  sliding  mode  func¬ 
tion  is  specially  designed  according  to  the  terminal  constraints, 
the  associated  sliding  mode  is  established  to  guarantee  the 
interception  with  the  specified  impact  angle  where  the  back- 
stepping  approach  is  utilized.  Moreover,  the  closed-loop 
dynamics  is  insensitive  to  uncertainties  and  disturbances  by 
the  virtue  of  the  strong  robustness  of  SMC.  The  guidance 
law  features  in  four  characteristics.  Firstly,  aerodynamic 


characteristics  are  taken  into  account  in  the  design  of  the  pro¬ 
posed  guidance  law.  As  a  result,  the  simulation  results  ob¬ 
tained  from  this  method  tally  well  with  the  actual  situation. 
Secondly,  there  is  a  wide  range  of  achievable  impact  angle  in 
the  usage  of  the  guidance  law.  It  guarantees  the  capability  of 
the  proposed  guidance  law  to  meet  different  mission  require¬ 
ments.  Thirdly,  the  guidance  scheme  is  quite  robust  to  initial 
conditions  as  well  as  parameter  perturbation.  Finally,  the  guid¬ 
ance  command  is  very  smooth  and  it  can  be  easily  realized  by 
the  attitude  control  system. 

2.  Problem  formulation 

In  this  paper,  a  two-dimensional  homing  scenario  is  consid¬ 
ered.  To  make  it  explicit,  the  engagement  geometry  in  inertial 
coordinate  frame  is  shown  in  Fig.  1.  The  missile  and  the  target 
are  represented  by  M  and  T;  the  missile’s  acceleration  is 
denoted  by  u ;  the  velocity  and  flight  path  angle  of  the  missile 
are  defined  as  V  and  y  respectively,  with  Vf  and  yf  being  their 
final  values.  In  addition,  the  missile  mass  is  represented  by  m 
and  the  gravity  acceleration  g  =  9.81  m/s  is  treated  as  a 
constant. 

The  equations  of  motion  in  the  terminal  guidance  phase  can 


be  given  as 

x  =  Vcosy 

(1) 

y  =  V  sin  y 

(2) 

D 

(3) 

V  = - g  sin  y 

m 

L  gcosy 

(4) 

7  “  mV  V 

where  the  position  coordinates  are  x  and  y.  The  terms  L  and  D 
are  the  aerodynamic  lift  and  drag  which  are  defined  by 

L  =  qSCL  (5) 

D  =  qSCD  (6) 

where  q  =  0.5pV2  is  the  dynamic  pressure  which  is  dependent 
on  the  air  density  p  and  the  velocity  V;  S  is  the  reference  area 
of  the  vehicle;  the  terms  CL  and  CD  are  the  lift  and  drag  coef¬ 
ficients  that  are  dependent  on  the  angle  of  attack  a  and  Mach 
number  Ma. 

Time  t  is  taken  as  the  independent  variable  in  the  equations 
of  motion  (l)-(4).  However,  the  flight  time  of  the  terminal 
guidance  phase  is  usually  not  the  key  question  we  are 
concerned  about.  In  Ref.  ,  the  downrange  x  is  used  as  the 
independent  variable  to  land  the  reusable  launch  vehicle 
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(RLV)  at  a  given  downrange  location.  Although  it  has  an 
advantage  in  dealing  with  such  landing  problems,  it  is  not 
appropriate  for  the  impact  angle  constrained  situations.  This 
is  because  each  independent  variable  only  possesses  one  corre¬ 
sponding  dependent  variable.  If  the  desired  impact  angle  is  set 
to  be  less  than  —90°,  each  x  will  have  two  corresponding  y 
when  the  value  of  x  is  larger  than  its  final  value.  As  a  conse¬ 
quence,  the  desired  impact  angles  which  are  less  than  —90° 
are  not  achievable.  Therefore,  a  new  independent  variable  is 
needed  to  handle  such  problems. 

The  perceived  wisdom  is  using  the  altitude  as  the  indepen¬ 
dent  variable  instead  of  the  downrange.  Because  the  altitude  of 
the  reentry  vehicle  is  usually  monotonically  decreasing  during 
the  terminal  phase,  each  altitude  has  only  one  corresponding 
downrange.  Note  that  the  independent  variable  used  in  the 
simulation  requires  a  monotonically  increasing  format.  There¬ 
fore,  a  new  variable  Y  =  y0  —  y  is  taken  as  the  independent 
variable  of  the  above  system  instead  of  t ,  where  y0  is  the  initial 
altitude  of  the  vehicle.  Thus,  the  system  model  can  be  changed 
to  the  following  equations: 


o 

o 

1 

II 

l-d 

II 

> 

(7) 

y  =  ^  =  -i 
-  dr 

(8) 

y  d  V  D  +  mg  sin  y 

d  Y  m  V  sin  y 

(9) 

,  dy  L  —  mg  cos  y 
^  d Y  mV2  sin  7 


The  time  becomes  an  incidental  variable  given  by 


d  Y  V  sin  y 


(10) 


(11) 


This  new  model  is  used  to  develop  the  guidance  law  detailed 
in  the  next  section,  and  the  derivative  of  the  flight  path  angle, 
/,  is  considered  as  the  control  input. 


3.  Impact  angle  guidance  law  design 


In  this  section,  an  SMC-based  guidance  law  is  designed,  which 
results  in  very  good  accuracy  both  in  miss  distance  and  impact 
angle  constraint. 

The  target’s  position  is  given  as  (xf ,  yf)  and  the  required  im¬ 
pact  angle  is  set  to  be  yf,  and  the  subscript  “f  ’  represents  the 
final  value  of  the  variable.  Since  yf  is  also  the  final  altitude 
of  the  vehicle,  we  have  Yf  =  y0  —  yf.  The  objective  is  to  enforce 
x  and  y  to  xf  and  yf  respectively  as  Y  reaches  Yf.  An  interesting 
result  can  be  obtained  from  Eq.  (7)  that  x'  only  relates  to  y, 
which  means  if  the  desired  impact  angle  is  given,  the  final  value 
of  x'  is  also  determined.  As  a  consequence,  the  impact  angle 
guidance  problem  can  be  solved  by  controlling  the  final  value 
of  x'. 

Motivated  by  this  idea,  the  sliding  function  is  designed  as 
follows: 

Si  =x-xf-x'(T-  Y{)  (12) 

Differentiating  Sj  with  respect  to  Y  yields 
s;=x'-x'  (13) 


It  is  obvious  that  at  the  time  of  interception,  the  value  of 
Y  —  Yf  will  go  to  zero.  Therefore,  if  Si  and  S\  also  go  to  zero 
at  the  same  time,  the  miss-distance  will  be  minimized  and  the 
impact  angle  constraint  will  be  satisfied.  To  make  the  first 
appearance  of  the  control  effort  /,  the  second  derivative  of 
Si  is  taken  and  one  gets 


S"  =  x"  =  — / 
sin  y 


(14) 


So  the  sliding  function  Si  has  a  relative  degree  of  two. 

In  Ref.  ’,  a  novel  second-order  sliding  mode  method  based 
on  an  adaptive  backstepping  approach  is  derived.  This  method 
can  be  used  to  ensure  that  both  S  and  S'  do  not  come  to  zero 
until  t  reaches  tf.  On  the  basis  of  this  idea,  a  similar  approach  is 
developed  which  is  able  to  guarantee  S  =  S'  =  0  to  come  into 
existence  at  Y  =  Yf  instead  of  t  =  tf. 


Theorem  1.  Consider  the  sliding  mode  function  S\.  If  S[  is 
treated  as  a  control  input  which  takes  the  form  of 


«>•  (15) 

both  S\  as  well  as  its  derivative  S[  will  go  to  zero  at  Y  =  Yf. 
Proof.  From  Eq.  (15),  one  gets 


,/  dSj nS i 

1  ~~dY~~  Yf-Y 


(16) 


Rearranging  variables,  Eq.  (16)  can  be  rewritten  as 

—  — n— -L— (jy  ,17) 


Assuming  the  initial  state  of  Eq.  (17)  is  (Tb,  Slb).  Then, 
integrating  Eq.  (17)  from  (Tb,  Sih)  to  the  later  point  (T,  Si), 
we  have 


P  Ps‘ = -  r 

Jslh  N  JYb 


Yf-Y 


-AY 


Thus,  the  following  equation  can  be  obtained 

1„ A  =  m  (*-*)’ 


(18) 


(19) 


Sib  "'{Y(-Yb)n 

Taking  the  exponential  of  both  sides  of  Eq.  (19),  the  analyt¬ 
ical  expression  for  Si  can  be  found  as 

Sib 


Si  = 


(Tf  -  Yby 


r  (Tf  —  Y)n 


(20) 


Differentiating  Eq.  (20)  with  respect  to  Y  yields 

From  Eqs.  (20)  and  (21),  we  can  conclude  that  Si  and  S\ 
will  both  go  to  zero  at  Y  =  Yf  as  long  as  the  value  of  n  is  larger 
than  one.  □ 

Therefore,  the  desired  trajectory  S\  is  defined  by  Eq.  (15), 
and  we  only  need  to  make  it  happen.  That  is  to  say,  a  control 
input  y'  needs  to  be  derived  to  enforce  S\  to  approach  and  then 
stay  on  the  desired  trajectory. 

A  new  sliding  function  is  defined  as 
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S2  =  g'+jt=Y’  n>l  (22) 

The  derivative  of  S2  can  be  obtained  as 


S'2  =  S"  + 


nS[  ( Ff  —  T)  +  mSj 

(jwf  : 


n  >  1 


(23) 


For  clarity,  let  M  =  1  ^  ^  n$l  we  have 

(Ff  —  Y)2 


^  =  -V/  +  M  (24) 

sin  y 

Next,  the  Lyapunov  theory  will  be  used  to  obtain  the  con¬ 
trol  input. 

Define  the  following  positive  definite  Lyapunov  function  as 

(25) 

Differentiating  V\  with  respect  to  Y  yields 


K  =  s2s'2  =  s2 


Xri  +  M 

sin  y 


The  control  input,  /,  takes  the  form  of 

y  7eq  7dis 


(26) 


(27) 


where  y'eq  is  the  equivalent  control  input  which  guarantees 
S'2  =  0  for  the  nominal  model,  and  y'dis  is  the  discontinuous 
switching  control  input  which  ensures  the  attractability  of 
the  sliding  surface. 

The  equivalent  control  input  is  selected  as 
?eq  =  —M  sin2  y  (28) 

Substituting  Eq.  (28)  into  Eq.  (27)  and  the  result  is  substi¬ 
tuted  into  Eq.  (26),  which  yields 

K  =  s2(^Yy'ih)  (29) 

\sin  y  J 

Then,  the  discontinuous  control  effort  is  chosen  as 
4is  =  —  ^  sin2  ysgn  S2,  k  >  0  (30) 

which,  on  substitution  in  Eq.  (29),  yields 


K  =  ~k\S2\  (31) 

If  S2  ^  0,  from  k  >  0,  one  has  Vx  <  0,  according  to  Lyapu¬ 
nov  stability  principle,  the  conclusion  that  the  system  is  asymp¬ 
totic  stable  can  be  drawn.  In  addition,  because  the  system  must 
perform  some  movement  along  the  trajectory  of  Eq.  (15)  to 
satisfy  the  terminal  constraints,  the  sliding  motion  should  be 
established  before  Y  =  Tf. 

Substituting  Eqs.  (28),  (30)  into  Eq.  (27)  and  the  result  is 
substituted  into  Eq.  (24),  one  gets 

Sf2  =  —ksgnS2  (32) 

From  Eq.  (32),  we  can  conclude  that  k  denotes  the  speed  of 
approaching  the  sliding  surface  S2  =  0.  Therefore,  the  switch¬ 
ing  gain  k  can  be  calculated  by 


£  =  iM01  ^  Yb  =pYf,  0<jp<l  (33) 

where  S2( 0)  is  the  initial  value  of  S2.  In  this  way,  the  sliding 
mode  S2  will  reach  zero  at  Y  =  Yh  and  keep  staying  on  the 
sliding  surface  thereafter. 


Till  now,  the  final  expression  of  the  control  effort  y'  can  be 
given  as 


7  = 


nS[  ( Tf  —  Y)  +  nS\ 

pwf 


\s2m 

pYf 


sgn  S2 


sin2  y 


(34) 


It  should  be  noted  that  the  guidance  process  is  over  at  the 
instance  the  value  of  Y  reaches  Tf.  In  other  words,  no  move¬ 
ment  is  allowed  along  the  trajectory  of  Sj  =  0.  Therefore,  un¬ 
der  the  circumstances  of  hitting  a  stationary  ground-based 
target,  the  available  range  of  yf  is  (— n,  0).  The  reason  of 
excluding  0  and  n  from  the  range  is  that  the  values  of  cot(0) 
and  cot(— n)  do  not  exist. 

Since  the  control  effort  will  enforce  the  system  states  on  the 
sliding  surface  S2  =  0  for  a  certain  period  of  time,  the  inherent 
chattering  problem  of  SMC  would  likely  to  be  experienced.  To 
solve  this  problem,  a  simple  saturation  function  sgmf(S)  = 
i  could  be  used  instead  of  the  sign  function,  where  s  is 

the  boundary  layer.  In  this  way,  the  system  is  no  longer  asymp¬ 
totic  but  could  only  converge  to  a  boundary  layer  about  the 
sliding  surface.  However,  by  setting  the  size  of  the  boundary 
layer  adequately  small,  the  error  can  be  limited  within  a 
negligible  range. 

In  many  scenarios,  the  desired  control  effort  for  a  lifting 
reentry  vehicle  is  the  angle  of  attack  a,  not  the  derivative 
of  the  flight  path  angle  y'.  However,  if  the  instantaneous 
value  of  y'  is  given,  the  angle  of  attack  can  be  obtained. 
Substituting  the  control  input  y'  into  Eq.  (10),  the  lift  L 
can  be  found  as 


L  =  mg  cos  y  —  m  V2y'  sin  y  (35) 

Then  from  Eq.  (5)  and  Eq.  (35),  the  lift  coefficient  CL  can 
then  be  obtained  as 

_  L  mg  cos  y  —  mV2y'  siny 

Ci=^= - ys -  (36) 

As  long  as  CL  is  achieved,  the  associated  angle  of  attack  can 
be  easily  derived  using  numerical  methods  such  as  Newton’s 
method  and  linear  interpolation. 


4.  Simulation  results 


In  order  to  verify  the  effectiveness  of  the  SMC-based  impact 
angle  guidance  law  advanced  in  the  previous  section,  numeri¬ 
cal  simulations  for  different  situations  are  presented  in  this 
section. 

In  the  simulations,  a  hypothetical  hypersonic  lifting  reentry 
vehicle  of  1000  kg  in  weight  and  2.5  m2  in  reference  area  is 
used  to  hit  a  stationary  ground-based  target.  The  aerodynamic 
model  of  the  vehicle  is  given  as  follows: 


Cl  —  CLo  +  CLad  (37) 

Co  =  CD0  +  KCl  (38) 

The  values  of  the  parameters  used  in  Eqs.  (37)  and  (38)  are 
given  below 

Clo  =  d—  arctan(10Ma  —  10)  —  0.035  (39) 

2071 

CLa  =  0.057e~°'654Ma  +  0.014  (40) 
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Table  1  Zero-lift  drag  coefficient  with  Ma. 


Ma 

Cdo 

Ma 

Cdo 

0 

0.0075 

5 

0.0092 

0.5 

0.0075 

5.5 

0.0087 

1 

0.0105 

6 

0.0082 

1.5 

0.0160 

6.5 

0.0080 

2 

0.0133 

7 

0.0078 

2.5 

0.0128 

7.5 

0.0079 

3 

0.0125 

8 

0.0080 

3.5 

0.0115 

8.5 

0.0080 

4 

0.0110 

9 

0.0080 

4.5 

0.0099 

9.5 

0.0080 

K=  1.85(1  _  e-°-2356M«)  (41) 

The  zero-lift  drag  coefficient  CD0  can  be  obtained  by  linear 
interpolation  among  the  data  in  Table  1  with  Mach  number  as 
the  independent  variable. 

The  atmospheric  properties  used  in  the  simulations  are 
based  on  the  1976  US  Standard  Atmosphere.  During  the 
implementation,  the  actual  control  effort  a  is  bounded  within 
(—30°,  +  30°).  First  of  all,  the  wide  capture  range  of  the  pro¬ 
posed  guidance  law  is  testified  by  the  simulations  with  different 
impact  angles.  Then,  some  cases  are  put  forward  on  the  basis 
of  vertical  impact.  Subsequently,  the  simulation  results  with 
fifth-order  system  delay  are  demonstrated.  Finally,  Monte 
Carlo  simulations  are  presented  to  validate  the  robustness  of 
the  guidance  law.  The  miss  distance  and  the  impact  angle  error 
are  defined  as  the  absolute  values  that  the  final  position  and 
terminal  flight  path  angle  deviate  from  their  desired  values. 

4.1 .  Simulations  with  different  impact  angles 

The  simulation  parameters  for  this  case  are  listed  in  Table  2. 
The  flight  trajectories,  flight  path  angle  histories  and  angle  of 
attack  profiles  for  this  case  are  plotted  in  Fig.  2(a)-(c),  respec¬ 
tively.  The  results  show  that  the  proposed  guidance  law  pos¬ 
sesses  a  wide  capture  range  of  impact  angles.  Fig.  2(c)  shows 
that  control  saturation  phenomenon  occurs  during  the  process 
of  achieving  the  design  goal  of  —150°  and  —170°.  This  is  be¬ 
cause  a  high  value  of  control  effort  is  necessary  to  ensure  the 
vehicle  to  maneuver  to  achieve  such  impact  angle  require¬ 
ments.  Besides,  the  miss  distances  as  well  as  impact  angle  er¬ 
rors  are  listed  in  Table  3.  The  miss  distances  are  all  limited 
within  1(T4  m  and  the  impact  angle  errors  are  no  bigger  than 


Table  2  Simulation  parameters  for  cases  of  different  impact 
angles. 


Parameter 

Value 

Parameter 

Value 

X0  (km) 

0 

Jt  (km) 

0 

Jo  (km) 

20 

Vf  (°) 

-10,  -30,  -90, 

-150,  -170 

V0  (m/s) 

1700 

n 

4 

To  (°) 

-10 

£ 

10“3 

Xf  (km) 

60 

P 

0.3 

- yf=-io°  - yf=-90°  — *-yf=-i70° 


0  20  40  60  80  100 

x  (km) 

(a)  Flight  trajectories  in  (*,>>)  space 


Fig.  2  Simulation  results  with  different  impact  angles. 


Table  3  Miss  distances  and  impact  angle  errors. 


Desired  impact  angle  (°)  Miss  distance  (m)  Impact  angle  error  (°) 


-10 

1.0262  x  10“5 

3.0327  x  10“6 

-30 

1.9299  x  10“5 

3.5163  x  10“6 

-90 

7.2760  x  10“u 

2.1587  x  10“n 

-150 

1.8934  x  10^5 

2.4075  x  10~6 

-170 

8.0306  x  10“6 

5.1661  x  10“6 

10  5(°).  The  results  show  excellent  performance  of  the  pro¬ 
posed  guidance  law  against  different  impact  angle  constraints. 

4.2.  Vertical  impact  cases 

In  this  part,  different  scenarios  of  vertical  impact  case  are  pre¬ 
sented.  In  all  these  simulations,  if  not  stated,  the  miss  distances 
are  limited  within  10-3  m  and  the  corresponding  impact  angle 
errors  are  less  than  10~4(°).  Additionally,  the  parameters  used 
in  the  following  cases,  unless  specified,  are  the  same  as  in 
Table  4. 

(1)  Case  A:  A  general  vertical  impact  case. 

In  this  case,  a  general  vertical  impact  case  is  presented  in 
Fig.  3(a)-(g).  The  vehicle’s  flight  trajectory  and  control  effort 
profile  are  plotted  in  Fig.  3  (a)  and  (d),  respectively.  For  most 
of  the  time,  the  angle  of  attack  changes  slowly.  However,  it  de¬ 
creases  quickly  during  the  final  period  of  the  flight  to  achieve 
the  desired  impact  angle.  As  the  flight  path  angle  is  nearly 
achieved,  it  will  converge  to  the  zero-lift  angle  of  attack.  As 
a  consequence,  the  altitude  of  the  vehicle  changes  gradually 
for  the  early  period  of  the  flight  but  it  descends  fast  when  it 
gets  close  to  the  target.  It  can  be  seen  in  Fig.  3(e)  that  the 
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Table  4  Simulation  parameters  for  vertical  impact  cases. 


Parameter 

Value 

Parameter 

Value 

x0  (km) 

0 

>Y  (km) 

0 

J’o  (km) 

20 

It  (°) 

-90 

Vo  (m/s) 

1700 

n 

4 

7o  O 

-3 

8 

10“3 

«o  (°) 

3 

P 

0.3 

Xf  (km) 

100 

0  20  40  60  80  100 

jt(km) 

(a)  Flight  trajectories  in  (x,y)  space 


(b)  Flight  path  angle  histories  (c)  Angle  of  attack  histories 


40  60 
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Fig.  3  Simulation  results  of  general  vertical  impact  case. 


sliding  mode  S2  occurs  at  the  instant  of  Y  =  0.3  Tf.  From  then 
on,  the  system  states  stay  on  the  trajectory  of  S2  =  0,  which 
guarantees  the  occurrence  of  Si  =  0  and  S[  =  0  at  the  final 
time  (see  Fig.  3(f)  and  (g)).  It  can  be  drawn  from  Fig.  3(b) 
and  (c)  that  the  impact  angle  requirement  is  satisfied  and  the 
velocity  history  is  satisfactory. 


Fig.  4  Simulation  results  of  vertical  impact  case  with  different  n. 


and  S\  as  well  as  the  flight  path  angle  histories  are  influenced 
apparently  by  the  variety  of  n.  It  is  because  the  behaviors  of  S\ 
and  S\  are  determined  by  the  value  of  n. 

From  Eqs.  (25)  and  (26),  we  can  see  that  the  value  of  n 
determines  the  shape  of  the  trajectories  of  S\  and  S\ .  For 
instance,  if  n  =  2  is  selected,  the  trajectories  of  S\  and  S\  are 
parabola  and  straight  line,  respectively.  Note,  from  Fig.  4(d) 
and  (e),  that  the  solid  lines  are  the  trajectories  of  S\  and  S\ 
when  n  is  set  to  2.  After  Y  =  Yh,  Eq.  (15)  is  satisfied,  and  thus 
the  two  trajectories  turn  out  to  be  parabola  and  straight  line 
respectively.  In  addition,  it  can  be  drawn  from  Eq.  (13)  that 
the  trajectory  of  S\  can  affect  the  behavior  of  the  flight  path 
angle  directly.  As  the  value  of  n  increases,  the  trajectory  of 
S[  tends  to  be  zero  more  quickly,  so  the  flight  path  angle  con¬ 
verges  to  its  desired  value  faster.  As  a  consequence,  less  time  is 
needed  for  the  vehicle  to  change  its  course  to  achieve  the  ter¬ 
minal  angle  requirement  (see  Fig.  4  (a)).  In  Fig.  4(c),  it  appears 
that  the  angle  of  attack  histories  are  smoother  and  their  peak 
values  become  larger  as  the  value  of  n  increases. 

(3)  Case  C:  Vertical  impact  with  different  p. 


(2)  Case  B:  Vertical  impact  with  different  n. 

In  this  case,  the  performance  of  the  proposed  guidance  law 
with  respect  to  different  n  is  tested  and  the  associated  results 
are  shown  in  Fig.  4(a)-(f).  The  values  of  n  in  the  simulations 
are  selected  as  2,  3  and  4,  respectively.  It  can  be  concluded 
from  Fig.  4  (b)  and  (f)  that  the  variety  of  n  does  not  have  much 
effect  on  the  velocity  histories.  However,  the  trajectories  of  S\ 


This  set  of  simulations  involves  implementing  the  guidance 
law  with  different  p.  The  simulation  results  are  described  in 
Fig.  5(a)-(f).  The  values  of  p  are  chosen  as  0.3,  0.5  and  0.7, 
respectively.  It  can  be  seen  from  Fig.  5(f)  that  the  sliding  mode 
occurs  at  the  desired  position  for  each  value  of  p.  From  the  sim¬ 
ulation  results,  it  is  obvious  that  varying  the  value  of  p  does  not 
affect  the  performance  of  the  proposed  guidance  law  apparently. 
Similar  to  the  previous  cases,  the  guidance  law  is  able  to  steer  the 
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Fig.  5  Simulation  results  of  vertical  impact  case  with  different  p. 


vehicle  to  attack  the  target  with  the  desired  terminal  angle  (see 
Fig.  5(a)  and  (b)).  It  can  be  noticed  from  Fig.  5(c)  that  a  larger 
magnitude  of  a  is  needed  for  a  larger  value  of  p. 

(4)  Case  D:  Simulations  with  fifth-order  system  delays. 


In  order  to  verify  the  performance  of  the  proposed  guidance 
law  in  practical  application,  the  simulations  in  this  part  are  car¬ 
ried  out  with  fifth-order  system  delays.  The  reason  for  choosing 
fifth-order  lag  is  to  make  a  clear  description  of  the  lag  model.  In 
this  canonic  model,  one  time  constant  represents  the  seeker, 
another  represents  the  noise  filter,  and  the  three  time  constants 
represent  the  flight-control  system.  The  fifth-lag  guidance 
system  possesses  the  following  transfer  function 


ar  1 

(t s  T  1 ) 


(42) 


where  t  is  the  guidance  system  time  constant,  and  ac  and  ar  are 
the  required  command  and  the  realized  command,  respec¬ 
tively.  The  time  constant  t  is  set  to  be  0.1s.  Other  parameters 
used  in  this  case  are  the  same  as  Table  2. 

Fig.  6(a)  and  (b)  show  that  the  vehicle  is  able  to  hit  the  target 
with  the  desired  impact  angle  even  under  high-order  time  delays. 
It  can  be  seen  from  Fig.  6(c)  that  the  time  histories  of  the  required 
command  and  the  realized  command  are  very  close  to  each 
other.  This  means  that  the  control  effort  obtained  from  the  pro¬ 
posed  guidance  law  can  be  easily  implemented  in  practical  appli¬ 
cations.  The  simulation  results  prove  the  robustness  of  the 
proposed  guidance  law  with  respect  to  high-order  system  delays. 


Fig.  6  Simulation  results  of  vertical  impact  case  with  fifth-order 
delays. 
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Fig.  7  Monte  Carlo  simulation  results. 


(5)  Case  E:  Monte  Carlo  simulations. 

In  the  previous  cases,  each  simulation  focuses  on  one 
situation.  In  this  case,  many  different  situations  as  well  as 
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high-order  system  delays  are  taken  into  consideration  in  one 
simulation.  Initial  perturbation  in  state  parameters  such  as 
altitude,  velocity  and  flight  path  angle  is  generated  randomly 
in  uniform  distribution.  The  range  of  dispersion  is  ±20%. 
Besides,  aerodynamic  modeling  error  is  applied  as  variation 
in  p  within  a  range  of  ±20%.  What’s  more,  random  noise 
ranging  between  —5%  and  5%  of  the  nominal  value  of  p  is 
added  to  the  aerodynamic  model.  The  Monte  Carlo  simulation 
for  500  runs  is  presented  in  Fig.  7.  From  Fig.  7(a),  we  can  see 
that  due  to  the  initial  perturbation  in  state  parameters,  the 
flight  trajectories  are  apparently  different  from  each  other. 
However,  all  of  them  converge  to  one  point  at  the  end  of  the 
flight.  Fig.  7(b)  and  (c)  show  that  the  results  are  equally  good 
even  under  such  severe  conditions.  For  majority  of  the  cases, 
the  miss  distances  and  the  impact  angle  errors  are  very  close 
to  zero.  The  miss  distances  in  all  cases  are  limited  in 
1.2  x  10-6  m,  while  the  impact  angle  errors  are  no  bigger  than 
5  x  10-5(°).  Good  robustness  of  the  proposed  guidance  law  can 
be  concluded  from  the  simulation  results. 

5.  Conclusion 

An  SMC-based  guidance  law  with  terminal  angle  constraint 
for  unpowered  lifting  reentry  vehicles  has  been  proposed  in 
this  paper.  A  special  sliding  mode  function  is  designed  on 
the  basis  of  terminal  constraints.  A  backstepping  approach  is 
used  to  obtain  the  control  effort  which  could  enforce  the  slid¬ 
ing  mode  function  together  with  its  derivative  to  zero  at  the 
time  of  interception.  The  guidance  approach  is  characterized 
by  its  simple  format  and  less  information  requirement.  The 
accuracy  of  the  guidance  law  has  been  demonstrated  in  the  pla¬ 
nar  engagement  scenario  with  considering  aerodynamic  char¬ 
acteristics.  A  wide  range  of  impact  angle  has  been  validated. 
In  addition,  different  simulation  parameters  are  adopted  to 
verify  the  robustness  of  the  proposed  guidance  law.  Moreover, 
the  performance  of  the  proposed  method  is  testified  by  imple¬ 
menting  the  guidance  law  with  respect  to  high-order  system  de¬ 
lays.  Finally,  strong  robustness  of  the  guidance  law  has  been 
drawn  from  the  Monte  Carlo  simulations.  The  future  work 
may  include  the  extension  to  three-dimensional  engagement 
scenarios  and  the  trajectory  shaping  ability. 
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